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Abstract 

A spherical tensor expansion of resonant inelastic x-ray scattering (RIXS) is carried out to express 
the RIXS cross-section in terms of a sum of products of fundamental RIXS spectra and tensors 
involving incident and scattered x-ray wavevectors and polarization vectors. Specific terms of this 
expansion are discussed in detail, in particular magnetic x-ray dichroism and powder samples. 

1 Introduction 

Resonant inelastic x-ray scattering (RIXS) is a remarkable tool to study electronic and phonon excitations 
in solids [T] . By selecting the suitable absorption edge one can choose a specific atomic species and orbital 
in a complex compound. The information content of RIXS is enhanced by the fact that the observed 
spectra depend on the polarization and direction of the incident and scattered beams. As a matter of 
fact, there is so much information in RIXS spectra that it is difficult to know whether a specific set of 
experiments exhausts it or not. The main purpose of our paper is precisely to determine a finite set of 
fundamental spectra in terms of which all possible experimental spectra can be expressed. More precisely, 
the RIXS spectrum obtained for a given wavevector and polarization vector of the incident beam (k and 
e) and scattered beam (k, and e^) is written as a sum of terms which are a fundamental spectrum [21 [3] 
(independent of the incident and scattered beam) multiplied by a given polynomial in k, e, k^ and eg. 

The number of available fundamental spectra decreases when the sample is a powder or when the 
polarization of the scattered beam is not measured. For example we show that the RIXS spectrum of 
a powder sample for electric dipole excitation and electric dipole emission with undetected scattered 
polarization is the sum of only three fundamental spectra. 

This approach was used with great success for the photoemission of localized magnetic systems |31 
[5l [6] and in x-ray absorption spectroscopy jH [3l [7] . Partial results were obtained in the fast collision 
approximation jSj or for the elastic scattering operator [SI (TUl [H] and general results were obtained in 
specific coordinate systems [T^l [131 E]. But, as far as we know, the fully general formula is new in 
coordinateless spherical tensor form. We showed for the case of x-ray absorption [7] that such geometric 
(i.e. coordinateless) expressions are useful (i) to disentangle the properties of the sample from those of the 
measurement; (ii) to determine specific experimental arrangements aiming at the observation of specific 
sample properties; (iii) to provide the good starting point for a group subduction leading to the angular 
dependence of specific crystal symmetries. 

Indeed, we show for example that magnetic dichroism might be observe on isotropic (i.e. without 
remanent magnetization) samples. 
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The range of validity of our final expression is rather broad. We only assume the validity of the 
Kramers-Heisenberg formula describing the scattering of light by a quantum system, where we neglect 
the Thomson scattering and the non-resonant term. We also assume that only the dipole and quadrupole 
transitions play a role, but this is an excellent approximation in the x-ray range [15j. Our result is a 
group theoretical consequence of the scattering cross-section. In particular, we do not assume that the 
states can be described by some LS or LSJ coupling, we do not make the fast collision approximation, 
we do not assume that some radial dipole and quadrupole matrix elements are independent of energy as 
in most other works on this subject. We even do not assume that a specific edge is measured. All these 
additional assumptions can be made in the final expression if necessary. 

We chose to give details of the computation to show that our result, even if it looks rather involved, 
can be obtained by straightforward angular-momentum coupling techniques. 

2 General case 

2.1 The Kramers-Heisenberg formula 

The scattering of light by a quantum system is described by the Kramers-Heisenberg formula, which was 
obtained before the advent of quantum theory [16] . The first quantum derivation of this formula is due 
to Dirac [H]. 

We start from the expression of the scattering cross-section: 

F N 

l^iFle- Pe-r|A.)(A.|.: . Pe— |/) ^^^^^ ^ _ _ ^^^^ 



m El — En — fv^s + ^7 

re is the classical electron radius: re — aoh/mc. \I), \N), \F) are respectively the initial, intermediate 
and final states. 7 is the total width of the intermediate state \N). The incident and scattered photons 
are characterized by the pulsation, wavevector and polarization vector w, k, e and a;^, k^, 63, respectively. 
Note that e* denotes the complex conjugate of es- 

In this expression, the first transition amplitude describes elastic scattering and will not be considered 
here. Ej being negative and large, Ei + huj can be small, and Ej — hujs is large. Hence, we can assume 
the third matrix element in Eq. [T] can be neglected with respect to the second one (although it could be 
treated with similar methods). Therefore, only the second amplitude transition remains in the expression 
of the scattering cross-section, yieding the well-known partial Kramers-Heisenberg formula: 

cTRixs = — — 7^ 7^ — — —■ S{Ef + fu^s- Ei -hw). 2) 

FN -1 iv I ; 

2.2 Multipole expansion 

First we transform the matrix element describing the absorption from the initial state |/) to the inter- 
mediate state I TV), (iV|e • Pe*''']/). 

To calculate {N\e ■ Pe''' ''|/), e'*' '' is expanded to first order: e'*' '' ~ 1 + ik • r. Hence, 

(iV|e-Pe^''-''|/) ^ (iV|e-P|/) +i(iV|k-re-P|/). 

The electric dipole matrix element is transformed by using the equation of motion of P which is P = 
{m/ih)[r,Ho]. l5^ Thus, 

(iV|e-P|/) = {m/ih){Ei-EN){N\e-r\I). 
For the quadrupole matrix element, one uses the identity from Ref. [18j . 

k-re-P = -(im/2n)[e • rk • r,i7o] + l/2(k X e) • L, 
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where L is the angular niomentum operator. One can neglect the angular momentum term (for an 
extended justificatfon, see Ref. [IS])- Therefore, 

(iV|6.Pe'''-"-|/) ^ --'-^{Ej~Er,){{N\e-r\I)+'-{N\e-rk-r\I)) 

1 

= -'-^{Ej-Er,)J2MN\e-r{k-rY\I), 

with /o = 1 and /i = |. 

Similarly, we transform the matrix element describing the emission from the intermediate state \N) 
to the final state (F|e* • Pe"'''=-''|7V): 

EF)(^{F\el-r\N}-^{F\e:-rks-r\N)) 

e'=o 

Finally, 



|iV) 



N 



im 
-^{Er, 



FN " ' i,i'=0 

xS{EF + hujs- Ei -huj). (3) 

In the following, we shall use: 



C 



E 



IFN 



Thus, Eq. [3] becomes: 



(El — En){En — Ep) 
Ei — En + huj + ' 
Se = S{Ep + hujs — Ej — huj). 



N f,f'=0 



Eq. |4]is a general expression describing the RIXS intensity for any combination of the absorption and 
emission transition operators. Each transition operator can be either pure electric dipole or pure 

electric quadrupole {E2), or a mixture fo both {Ei + £^2)- 



2.3 The Kramers-Heisenberg formula expressed in spherical tensors 

In this subsection, the expression of the RIXS cross-section (Eq. |3]) is transformed using spherical tensors 
and their coupling properties. For a short introduction to spherical tensors and their application to x-ray 
spectroscopics, the reader is referred to Ref. [7] for the case of the X-ray absorption cross-section. Let 
us first briefiy explain the notation. A rank spherical tensor T is written as T^^^. This must not be 
confused from with is the ^-th power of T. 

Cartesian vectors, such as e, r or k are written in their usual form, i.e., without brackets. However 
one shall keep in mind that Cartesian vectors correspond to first-rank spherical tensors, and as such they 
shall indifferently be written as e^^\ r'^^\ k^^^. 
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2.3.1 Expression of e • r(k • r)^ 

Here we show that e • r(k • r)^ can be expressed as ge< {e (X) k^}(^+i) I where gg is a constant. 

For £ = 0, according to Ref. [7 (Eq. 13): 

e-r = -V3(--^eT) = -\/3{e®r}0 = -\/3|{e®k"}(i)®r(i)}*°\ 
For i = 1, according to Ref. [7 (pp.5): 

e-rk-r = V5{{e k}(2) ,g, {r « r}(2)}(") ^ ^/5|{e (g) k^j'^) ® r^^) j^"^ 

Thus, e-r(k-r)^ = c/^|{e k^}(^+i) r(^+i)|^°\ with r(i) = r, go = -^3, r^^) = {r(8)r}(2) and 

51 = Vs. 

Then, Eq|4] becomes: 

fRIXS 

(0), , , ,r, , ^ 



F N £,£'=0 

^'^^E |E^^^^^ E Mi'mi'{{e^^T+'^ « (7V|r(^+i)|/)}^°^{K «kf ® (F|r(^'+i)|iV)}^°^ 

F N e,e'=o 



F ' N e,e'=o 
with hi = /^5£: /iq = -a/3, /ii = f a/5. 



If one defines rj^|^^ = (iV|r(^+i)|/) and r^^^^^ = (F|r(^'+i) |iV), one gets: 



2.3.2 Recoupling spherical tensors 

Here we use the recouphng identity (Ref.fT, Eq.l4), which is: 

ff=|a-<i| V(2a+ l)(2rf+ 1) 

where the scalar product • of two spherical tensors P'^-' andQ^^' is defined by P^^^-Q^^^ = J2^=-gi~'^V P^^^ Q~\- 
Eq. |6] becomes: 

^ E I E E 3) « « " « ^''""'}'" ■ " « 

where g runs from |^ — to (i + i' + 2), £ and ^' run from to 1. 
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The first tensor product |{e* (8)kf ® {e(g)k^}'^+^) | gathers the characteristics of the incident 

beam (e, k) and of the scattered beam (es,ks). The variables dealing with the sample are gathered in 

the second tensor product Ir^/^^' ® r^^^'j 



(9) 



Then, defining 



N 



_ ^ ^ (-^/ ~ Em){En - Ep) ,je'+i) ^ ji+i) 

N 

one gets 



F g,' 



2.3.3 Expansion of the square 

Here we square the amplitude inside the sum over the final states in Eq. 1101 
When expanding the square, Eq. 1101 becomes: 



crixs = ^Oe / / — , 

F V(2^i + 3)(2f: + 3)(2£2 + 3)(24 + 3) 

(51 



{{6, ® klH(^^+i) ® {e* ® k^^}(^^+i)}^''' . 4|)(£2,4)*, (11) 
where Qi runs from \ti — 1[\ to {ti + £[ + 2), with li and l[ equal to or 1. 



In the Appendix, we show that 4|\4,4)* = ^^^^^(^2,4). 
Therefore, 

crixs — C5e/ / — , ^ 

v/(24 + 3)(2fi+3)(2^2 + 3)(2^^ + 3) 



{{e,0k?}(^^+i)0K0k^^}(^^+i)}^''^4^;\£2,4) (12) 



If one defines 



- {{e:®k5}W+i)0{e0k^^}(^^+i)} , (13) 

X^^^) = {{e,®k'i}(^i+i)®{e*®k^^}(^^+i)}^''\ (14) 

n,,+i,,,+i,,i+i,,^+i = ^(2^1+3) (2fi + 3) (2^2 + 3) (2f2 + 3) , (15) 
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one gets 



p » », -l-l£i + l,£2 + l.f'i + l,i^2 + l 



31 1 

32, ^2, ^'2 



Using the identity P^'^' • Q^"^ = {-lY^2a + l {P(°) (g) Q(°)}(o) (pp.64-65 of Ref. [Tg'), we have: 
and then: 

^ ^ {~lY'+''^-'->^h,,hl,hlh,^ ^(251 + 1) (232 + 1) 

CTRIXS = C()£;2^ 2^ 

ff2!^2!^'2 

{xj^^^ ® 4^1^(4, ri)}(°H^2''^ ®4'j^\^2, 4)}*°'- (19) 

Using Eq. [71 we have: 

{x}^^) ®4])(^i,£'i)}(o){x(3^) ®4f^(^2,4)}'"^ 

91+92 ( —^'\'^ 

= 5] -_44=_{x^)®x(^^)}('^).{4|)(^i,fi)®4^;\£2,4)}('^) (20) 

a=\gi-g.\ V (251 + 1)(2.92 + 1) 

Then, 
ctrixs 

/ 2^\i?2+^2^P2 ^ ^ 

= c5eY. —r ^^^iLA ^ (-i)nx^)«x(^^)}W.{4|)(£i,f,)«4^;)(£2,4)}(''\ 

F g,,l^/^ ll£l + l,£2 + l,<?i + l,<?i + l 

92,^2,^2 

with a running from |gi — 32! to gi + (72- 
In order to simphfy the tensor product: 

one uses the following identity (pp.70 of Ref. [19 ): 

(fe) fa 6 c"! 

where \a — d\ < g < a -\- \b — e\<h<b-\-e and \g — h\ < k < g + h. 
This yields: 



7; + l ^1 + 1 gi 

{x^)®x(^^)}W = 5]ng,,g,,5.J4 + i ^2 + 1 52^ 
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c a 



with - {{{e:®kfi}«+i)®{e,®kf^}(^2+i)}(^)0{{e,^k^i}(^i+i)0{e*,^k^2}(f2+i)}(-)|*"\ 

According to the triangular conditions, the 9j-factor is zero if any of the fohowing conditions is not 
satisfied 

\b-c\ <a< b + c . 

2.3.4 Simplification of {X}^'') 

In {X}^'^\ the variables concerning the scattered beam (es,ks) are gathered in the first tensor product, 
while those concerning the incident beam are gathered in the second product. Therefore in order to 

simplify {X}^'^), one shaU write it as {X}^"') = (g) {Xabs^'^'^y''^ and then simplify {X^^}^''') 

and {Xabs}'''^^- According to Eq. [5T1 one gets: 



fl 4 ^1 + 1] ^ 



U V c 



'I ll ll + l 



U V 



(&) 



U V 



'I ^1 + 



The 9j-factors vanish if any of the following triangular conditions is not fullfilled: 

<u< 2, 
1^1-41 <v< £1+4, 
< u' < 2, 

|£'i-4l <^^'< ^'1+^2, 

2.4 Final expression of ctrixs 
Finally, one obtains: 

Ctrixs 



n ^1 ^i + n n £'i + n r^'i + i ^1 + 1 gl^ 

F 91,92 a,b,c ^ ^ U V c ) ^u' v' b ) ^ b c a) 



^i'^2 v,v' 
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with 



■0 



-V3 



C 




{aY'''> 



ft2 UJ 

{{e (g) e*}(") ® {k^i (g) k^^ 
{{e:®ej("')®{kfi®k'/}("')} 




Af >{£,£') 





(Ar|r(^+i)|/), 
(F|r(^'+i)|Ar). 



In Eq. [22] the RIXS cross-section is written as a sum of fundamental spectra, < Apj {ii,£[) <Si 



Ajp (^2,^2) f I which are only sample-dependent for a given couple of transition operators. The fun- 



damental spectra are weighted by coefhcients which depend only on the experimental conditions, i.e., 
transition operator, type of polarization and experimental geometry, in the absorption as well as in the 



emission. The |{£;}W(8){A}(^)1 weights are polynomial in e, e.,, k and kg. The three 9j-factors (and 



the related triangular conditions) ensure that these variables are coupled in the RIXS process in a correct 
(physical) way. 

Applications of the general expression to particular cases will be given as examples in the next section. 

2.5 Summary: Possible values of the different indices 

Here we summarize the possible values taken by the different indices in Eq. 1221 

9 < £i < 1 (for i = 1,2 ): these indices deal with the transition operator in the absorption. 
£i = for dipole excitation, £i = I for quadrupole excitation 

• < I'i < I (for i — 1,2 ): these indices deal with the transition operator in the emission. = 
for dipole emission, £^ — 1 for quadrupole emission 

• \£i — £'^\ < gi < £i + f.^ + 2 { ior i = 1,2 ) : these indices couple absorption and emission transition 
operators. 

• < u < 2: u deals with the polarization of the incident beam. 

• < u' < 2: : u' deals with the polarization of the scattered beam. 

• \£i — £2] < V < £1 + £2- V deals with the direction of the incident beam. 

• \u — v\<c<u + v and \£i — £2] ^ c < £1 + £2 + 2: c gathers all characteristics of the incident beam. 

• ~ ^21 — — + ^2- deals with the direction of the scattered beam. 

• \u' — v'\ < b < u' + v' and ]£[ — £'2] < b < £[ + £'2 + 2: b gathers all characteristics of the outgoing 
beam. 

• \b — c\ < a < b + c and I51 — 52I < a < .9i + 52 : a couples everything. 
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3 Particular cases 



In this section we consider several particular cases. 

3.1 Conditions due to the type of polcirization 

When the incident beam is circularly polarized, m = 0, 1, 2. When the incident beam is linearly polarized 

light, fi = or 2. Note that there is a basic difference between the polarization of the incident beam 
and that of the scattered beam. On the one hand, the incident beam is prepared by a synchrotron 
radiation source and monochromators. Therefore, e can be tuned as one wishes. On the other hand, the 
polarization state of the scattered beam is entirely determined by the incident beam and the sample. It 
cannot be tuned. However, we can measure some polarization properties of the scattered beam. If we 
want to measure the intensity along the polarization direction e,,, then we introduce Cg in the formula. If 
we do not measure the polarization state, then we take the trace of the density matrix representing the 
polarization state of the scattered beam. In other words, the measured cross-section is equal to the the 
sum over two perpendicular directions of the scattered polarization (see Appendix 1). In this case u' = 1 
is zero, since ({e e*}^^^) = 0. 

3.2 Electric dipole excitation, electric dipole emission 

We consider the case of an electric dipole transition in the absorption followed by an electric dipole 
transition in the emission : ii = 0, £2 = 0, i\ = and £'2 = 0. Thus, we have: 

< 51 < 2, 
< 32 < 2, 
< a < 4, 
0< & < 2, 
< c < 2. 

Since £1 = 0, £2 = 0, we have v = 0. This implies c = u. Since £\ = and £'2 = 0, we have v' = and 
b = u' . These conditions allow to calculate the values of some 9j symbols. 







h + l] 














Pi 




V 






[cOcj 







^i + l] 














Pi 




v' 






[bObj 



3V25+T' 



Thus, we have: 

fi 1 rn 

F 91,92 a,b,c [ 6 C a 



Additionnally, 



{eY'^ = {{e:®ejW®{kO®kO}(o)}^''^ = {e>e,}W 
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3.3 Electric dipole excitation and emission, powder sample 

A further simplification arises when the sample is a powder. The ensuing isotropy implies that a = 0. 
Thus, we can use 



5/0/ v'(25i + l)(26 + l) 111 b 

This implies that b = c and gi = 92- Therefore, 



gfi ^ ^ 

{A^fiM ^A^fMY^ 5{Ef + ^s-Ei-nw) (24) 



• For the term 6 = 0: 



with g = 0,1,2. 



Thus, 

aRixs(6 = 0) = V^|±l ^ {4s](0,0) ® 4F(0,0)}^°^.5(i?f + huj, - Ej - huj). 

3=0 F 



For the term 6=1: 



ml _ fll 2l 1 
lllf " llllj 6' 



11 01 ^ _1 
lllf 3' 



{e®e*}« = Aexe* = -^k, 
\/2 \/2 

where Pc is the rate of circular polarization. Note that Pc is positive for a right circular polarization 
in the traditional sense (i.e. for a negative helicity). Pc = if the polarization is linear. 

= -l=P,P,,,k,.k=^(|e.e:p-|e.e.n. 

This result is interesting. We see that, for an unoriented sample, in particular for a sample without 
permanent magnetization direction, some circular dichroism can be observed if we measure the 



10 



circular polarization of the scattered beam. However, we need to avoid the configuration where 
the incident and scattered x-ray wavevectors are perpendicular. This remark shows the power of 
the geometric (coordinateless) approach. Note that, in this expression Pc.s does not denote the 
degree of circular polarization of the scattered beam, but the efficiency of detection of the circularly 
polarized beam. Note also that e • e* and e • are the only non-trivial scalars that can be built from 
e, Cg and their conjugates. Therefore, we could expect them to be involved in this formula. 

• For the term b = 2: 

(11 g\ ^ 4 
1112/ (2-g)!(3 + 5)!' 

with g = 0,1, 2. We have also 

If the polarization of the scattered beam is not detected, then the corresponding term 6 = 2 is 
calculated using the relation: ({eg (g) ej}*-^^) = — ki^V^ (see Appendix). Thus, 



3.4 Electric quadrupole excitation, electric dipole emission 

Here we consider the case of an electric quadrupole transition in the absorption followed by an electric 
dipole transition in the emission: = 1, £2 = 1, ^'i = and £3 = 0. Thus, we have: 



< 51 < 3, 
< 52 < 3, 
0< a< 6, 
< 6 < 2, 
< c < 4. 

Since £1 = 1 and £2 = 1, we have 0<v<2. Additionally, w ^ 1 since {k^ (g) k^jfi) = -^k x k = 0. Thus, 
= or 2. 

Since £[ = and £'2 = 0, we have v' = and b = u' . 



£1 


tl + l] 






£2 








V 

















(-'2 














\bOb] 



The second 9j-factor is: ^ 1 ^2 

^v! v' 

7i + l £1 + 1 51 
• The third 9j-factor is: ^ ^2 + 1 ^2 + 1 92 

b c a 



3.4.1 Powder sample 

Isotropy implies that a = 0. This implies that gi = 52 and b = c = u'. 









1252 




[ & c a j 
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3.4.2 Circularly polarized light 

For circularly polarized light we have 



1 1 2' 

The first Qj-factor is: ^112 

Ival 



4 Conclusion 

We have presented a detailed derivation of the RIXS cross section in geometric spherical tensor form. 
Some special cases were studied to illustrate the ability of this expression to describe global properties of 
the sample (such as being powder). 

We have also clarified the meaning of the scattered beam polarization vector as a way to interpret a 
measurement of that beam. 

This general formula is ideally suited for a group subduction corresponding to particular symmetry 
groups. This will be illustrated in a forthcoming publication. 



5 Appendix 

One of the most powerful features of geometric expressions is that they can be calculated in a specific 
coordinate system adapted to a particular problem, and then be valid in any system. We illustrate this 
by describing the average of the coupling of e and e* : we show that 



({e«)e*}W) 



1 

= 0, 

We take a reference frame where k is along Oz. The linear polarization vector is e = {cos ip, sin ip,0) and 
the corresponding spherical tensor components are = c'q^ = 0. Therefore, it is easy to 

calculate 



(2) _ 

2 ' 

(2) 



{e(3e*}'^> = 0, 



(2) 
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The average over two perpendicular polarizations {ip and tp + 7r/2) or the average over all ip gives us 

{{e^e*}^il) = 0, 
{{e^e*}^il) = 0, 

In that reference frame wc also have knJ = S,n,Q^/V3. Thus, ({e® e*}^^)) = -k(2)/2. Moreover, this 
relation is true in all reference frames because it is a relation between two tensors. 
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We can also consider elliptically-polarized x-rays. Thie most general polarization vector in a frame 
where k is along Oz is 

cos X cos V' + i sin X sin ij) 
e = I cos X sin ^ — i sin x cos i/' 


for which the degree of circular polarization is sin 2x- In particular, x = 7r/4 for a fully circularly polarized 
x-ray. Then, 

{e®e*}« = 0, 
{e®e*}« = --Lsin2x, 

and 







= ^cos2xe^ 






= 0, 






1 



Therefore, by using left and right fully circularly polarized beams, we obtain the same average as with 
linear polarization. 



6 Complex conjugate 



We want to calculate the complex conjugate of The operators x, y and z are Hermitian. 

Therefore 

(iV|a;|7) = (iV|a;t|7) = (7|a;|iV)*, 

and (A^|a;|/)* = (/|x|A^). The same is true for y and z. For the corresponing spherical tensors, we have 

(7V|r«|/>* = -(l/v^)((^|x|/)+i(7V|t/|/))* = -(l/V2)((J|a;|iV)-i(/|y|7V)) =-(7|rW|iV)*. 

An analogous calculation for the other components of r^^^ gives us (ATlr^^^ |7)* = {—l)^{I\r^_^}^\N). For 
£= 2, we have 

AA' AA' K 

AA' K 

The Clebsch-Gordan coefficient implies that A + A' = /i. We replace A and A' by — A' and — A, respectively, 
and we obtain 



(Ar|rf|7)* = {-irJ2(^-X'l-X\2f,)J2{I\r'-x^\K){K\r^,^\N) 

AA' K 

= {-ir -^'1 -A|2/.)(7|rWrl^) \N}. 



AA' K 
AA' 

Now, we use the symmetry (^imi^2TO2K3W3) = {£2 — ^2^1 — "t-iKs — "t-s) to obtain 



(7v|r(2)|7)* = (-irE(i^i^'i2-/^)(^i^i'^^i'^i^) = (-ir(^k->)- 



AA' 

13 



A recursive use of this argument leads to {N\rm\I)* = (— l)^(/|r^^|A^) for any £. 
We proceed with 

m'm 

The complex conjugate of is 

X; = Y.{i'm'£m\gj){F\r^^}\Nr{N\rii^\ir = ^{£'m'£m\gj){-ir+"^\N\/!l\F}{I\r'^l\N}^ 

m'm m'm 

(2.) 

The same reasoning as for tlic complex conjugate of {N\rli \I) gives us 
Finally, 

a(9) (p fiy _ / V ~ -^^)(^^ ~ rr(^+i) e- r^^'+^)^(a) 



^ Ei - En + fvjJ - 



(25) 



Therefore, if we define 



^ El - En + huj - 



we have )* = f ) 

In applications, we compute only {N\r^^^\I) and not (/jr^^'' jiV) that we need to calculate Ap] ^{£,£')*- 
However, we can easily obtain it by (7V|r^^^|/) = (-l)^(/|ri^]^|iV)*. 

Wc also have to calculate the complex conjugate of |{e*(g)kf +^)(g){e{X'k^}('^+^)| . The calculation 

is very similar. Consider a vector z = a + ib, where a = {ax,ay,az) and b = {bx,by,bz) are real. We 
have 

(z^^^)* = --^{ux + ibx + i{ay + iby))* = --^{ttx - ibx - itty - by)) 

= --^(a^-ifex-iK-iM) =-(z*)-i- 

The calculation of the other components gives us (z^^^)* = (— )'^(z*)^|. Therefore, (ex)* = (— l)'^(e*)_A 
and {{€*)x)* ^ (-l)^(es)_A. The vector k is real and (k^')* = (-l)'"kL^^. By repeating the proof for 

A^^l{£,£') we have 

({e0k^}(^+i))* = {-l)r-J2i^-Xi-i,\i+lm)elkW. 

Since we do not interchange lA and £ii, we must use the symmetry relation 

(4mi4m2|4m3) = -mi4 -m2|4 -mg), (27) 

to get 



({e®k^}(^+i))* = i-l)"^^{lX£i,\£+l-m)elkl'^ = i-l)"\e*®k'}':!i'\ 

Xfj. 
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Here, the symmetry relation docs not bring any additional sign because £3 — ii + ^2- Similarly, |^{e* (8) 
kf }^'+^^)* = (-!)'"'{£« i^U-'r^'^^- We now do the same calculation for {{e* (g) kf (g) {e ® 

k^}(^+i) I ^ but here, the sign (—1)^+^ f must be retained and we get 

({{e: ® kf ® {e ® k^}(^+i)}^'y = (-l)7(_i/+^'-fl{{e^ ^ kf ® {e* ® k^}(^+i)}|_'\ 

Finally, the complex conjugate of X = |{e* (g kf O {e (g k^}(^+i)|^^^ • ) is 

X* = (-l)^+^'-s{{e, ® kf {e* ® k^}(^+i)}^'^ . 4f (^,/). 
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